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Abstract 

The charged particle production in proton-proton colhsions is studied with the LHCb 
detector at a centre-of-mass energy of y/s = 7TeV in different intervals of pseudorapid- 
ity 7]. The charged particles are reconstructed close to the interaction region in the vertex 
detector, which provides high reconstruction efficiency in the rj ranges —2.5 < t] < —2.0 
and 2.0 < rj < 4.5. The data were taken with a minimum bias trigger, only requiring one 
or more reconstructed tracks in the vertex detector. By selecting an event sample with at 
least one track with a transverse momentum greater than 1 GeV/ c a hard QCD subsample 
is investigated. Several event generators are compared with the data; none are able to 
describe fully the multiplicity distributions or the charged particle density distribution as 
a function of rj. In general, the models underestimate the charged particle production. 



Keywords: minimum bias, underlying event, particle multiplicities, LHC, LHCb 



^Authors are listed on the following pages. 



ii 



LHCb Collaboration 



R. Aaij23, C. Abellan Beteta^^'", B. Adeva^^, M. Adinolfi^^^ C. Adrover^, A. AfFolder^^ 

Z. Ajaltouni^, J. Albrecht^'', F. Alessio^'^, M. Alexander'''', G. Alkhazov^^, 

P. Alvarez Cartelle^^, A. A. Alves Jr^^, S. Amato^, Y. Amhis^^, J. Anderson^^, R.B. Appleby^°, 

O. Aquines Gutierrez-^", F. Archilli^'^'^'^, L. Arrabito^^, A. Artamonov M. Artuso^^'^'^, 

E. Aslanides^ G. Auriemma^^'"^, S. Baclimann", J.J. Back^^ D.S. Bailey^^, V. Balagura30.37^ 

W. Baldinii^, R.J. Barlow^o, C. Barschel^'^, S. BarsuF, W. Barter^^^ Bates'''^, C. Bauer^o, 

Th. Bauer23, A. Bay3^ I. Bediaga^ S. Belogurov^o, K. Belous^^, I. Belyaev^o.s? ^ 

E. Ben-Haim^, M. Benayoun^, G. Bencivenni^^ , S. Benson^^, J. Benton^^, R. Bernet^^, 
M.-O. Bettler^^, M. van Beuzekom^^, A. Bieni\ S. Bifani^^ T. Bird^°, A. Bizzeti^^'^, 
P.M. Bj0rnstad5O, T. Blake^^ F. Blanc^s, C. Blanks''^ J. Blouwi\ S. Blusk^^^ A. Bobrov^^, 
V. Bocci^^, A. Bondar^^, N. Bondar^^, W. Bonivento^^, S. Borglii'^'^'^^, A. Borgia^^, 

T.J.V. Bowcock'*^, C. Bozzi^^, T. Brambach^, J. van den Brand^^, J. Bressieux'^^ , D. Brett^'^, 

M. Britsch^°, T. Britton^^, N.H. Brook''^, H. Brown^^, A. Biichler-Germann^^, I. Burducea^^, 

A. Bursche^s, J. Buytaert^^, S. Cadeddu^^ O. Callot^ M. Calvi^Oj, M. Calvo Gomez^^'", 

A. Camboni^^, P. Campana^^'^'^, A. Carbone^'', G. Carboni^^'*^, R. Cardinale^^'*'^'', 

A. Cardini^^, L. Carson^^, K. Carvallio Akiba^, G. Casse^^, M. Cattaneo^^, Ch. Cauet^, 

M. Charles^i, Ph. Charpentier^^, N. Chiapolini^s, K. Ciba^^, X. Cid VidaP^, G. Ciezarek'^^, 

P.E.L. Clarke'"5.37^ M. Clemencic^^ H.V. Cliff^^^ J. Closier^^, C. Coca^^, V. Coco^^, J. Cogan^ 

P. Collins^^, A. Comerma-Montells^s, F. Constantin^s, G. Conti^^, A. Contu^\ A. Cook''^, 

M. Coombes^^, G. Corti^^, G.A. Cowan^^, R. Currie""^, B. D'Almagne'^, C. D'Ambrosio^'^, 

P. David^ P.N.Y. David^^, I. De Bonis^ S. De Capua^^'*^, M. De Giants, F. De Lorenzi^^, 

J.M. De Miranda\ L. De Paula^, P. De Simone^^, D. Decamp^, M. DeckenhoflP, 

H. Degaudenzi^^'^'', M. Deissenroth^^, L. Del Buono^, C. Deplano^^, D. Derkach^'^'^'^ , 

O. Deschamps^, F. Dettori^^, J. Dickens^^, H. Dijkstra'^^, P. Diniz Batista^, 

F. Domingo Bonal^^'", S. Donleavy^^, F. Dordei^^, A. Dosil Suarez^^^ D. Dossett^^, 
A. Dovbnya''°, F. Dupertuis3^ R. Dzhelyadin^^, A. Dziurda^^, S. Easo^^ U. Egede^^, 
V. Egorychev^o, S. Eidelman^^, D. van Eijk^^, F. Eiselei\ S. Eisenllardt''^ R. Ekelhof^, 
L. Eklund^^, Ch. Elsasser^^, D. Elsby^^, D. Esperante Pereira^^, L. Esteve^^^ 

A. Falabellal6'l^'^ E. Fanchini^O'J , C. Farber", G. Fardell^^ C. Farinelh^^, S. Farry^^ 

V. Fave^^, V. Fernandez Albor^^, M. Ferro-Luzzi^'^ , S. Filippov^^, C. Fitzpatrick^^ , 

M. Fontana^°, F. Fontanelh^^'*, R. Forty^^, M. Frank^^, C. FreF, M. Frosini^^'-^'^T^ Furcas^", 

A. Gallas Torreira^^, D. Galh^^''^ M. Gandelman^, P. Gandmi^i, Y. Gao^, J-C. Garnier^^ 

J. Garofoh^^, J. Garra Tico^^^ Garrido^^, D. Gascon^^, C. Gaspar^''', N. Gauvin^^, 

M. Gersabeck37, T. Gershon'^^'^^, Ph. Ghez^ V. Gibson^^ V.V. Ghgorov^^, C. Gobel^^, 

D. Golubkov^o, A. Golutvin'^^'^O'Sr^ a. Gomes^, H. Gordon^i, M. Grabalosa Gandara3^ 
R. Graciani Diaz^^, L.A. Granado Cardoso^'', E. Grauges'^^, G. Graziani^'', A. Grecu^^, 

E. Greening5\ S. Gregson^^^ B. Gui^^^ E. Gushchin^^, Yu. Guz^'', T. Gys^^ G. Haefeli3^ 

C. Haen^^, S.C. Haines^^, T. Hampson'^^, S. Hansmann-Menzemer^^, R. Harji^^, N. Harnew^^, 

J. Harrison^'^, P.F. Harrison^^, J. He^, V. Heijne^^, K. Hennessy^^, P. Henrard^, 

J. A. Hernando Morata^^, E. van Herwijnen^^, E. Hicks^^, K. Holubyev^^, P. Hopchev'^, 

W. Hulsbergen23, p. Hunt^i, T. Huse^^ R.S. Hustoni^, D. Hutchcroft^s^ D. Hynds^^ 

V. lakovenko^^, P. Ilten^^, J. Imong^^, R. Jacobsson'^^, A. Jaeger^^, M. Jahjah Hussein^, 

E. Jans^3, F. Jansen^^, P. Jaton^^, B. Jean-Marie^, F. Jing^, M. John^^, D. Johnson^^, 



iii 



C.R. Jones'^^^ B. Jost^^, M. Kaballo^ S. Kandybei^°, M. Karacson^^, T.M. Karbach^, 

J. Keaveneyi^^ I.R. Kenyon5^ U. Kerzel^^ T. KeteP, A. Keune^^, B. Khanji^ Y.M. Kim^^ 

M. Knecht^^, P. Koppenburg^^, A. Kozlinskiy^^ , L. Kravchuk^^, K. Kreplin^^, M. Kreps^^, 

G. Krocker^^, P. Krokovny^^, F. Kruse^, K. Kruzelecki^'', M. Kucharczyk^"'^^'^ '''■^ , 

T. Kvaratskheliya^O'^^, V.N. La Thi^^, D. Lacarrere^^, G. Lafferty^^, A. Lai^^ D. Lambert^^ 

R.W. Lambert^^, E. Lanciotti^^, G. Lanfranchi^^ , C. Langenbruch^^, T. Latham^'^, 

C. Lazzeroni^^, R. Le Gac^, J. van Leerdam^^, J. -P. Lees^, R. Lefevre^, A. Leflat^^''^^, 

J. Lefrangois^, O. Leroy^, T. Lesiak^^, L. Li^, L. Li Gioi^ M. Lieng^ M. Liles^^ R. Lindner^^, 

C. Linn^^ B. Liu^, G. Liu^^, J.H. Lopes^, E. Lopez Asamar^^, N. Lopez-March^^ , H. Lu^^'^, 

J. Luisier^^, A. Mac Raighne^^, F. Machefert^, LV. Machikhiliyaii^'^° , F. Maciuc^°, 

O. Maev29'37, j. Magnm\ S. Malde^\ R.M.D. Mamunur^"^, G. Manca^^-'^, G. Mancinelli'' , 

N. Mangiafave^^, U. Marconi^'', R. Marki^^, J. Marks^\ G. Martellotti^^ , A. Martens^, 

L. Martin^\ A. Martm Sanchez'^, D. Martinez Santos^'^, A. Massafferri\ Z. Mathe^^, 

C. Matteuzzi^o, M. Matveev^^, E. Maurice'^, B. Maynard^^^ a. Mazurov^^'^^.s?^ 

G. McGregor^", R. McNulty^^^ C. Mclean^"^, M. Meissneri\ M. Merk^^, J. Merkel^ 

R. Messi^i'*^, S. Miglioranzi^^ , D.A. Milanes^^'^'^, M.-N. Minard'', J. Molina Rodriguez^'', 

S. Monteil^, D. Moran^^^ Morawski^^, R. Mountain^^, L Mous^^, F. Muheim^*^, K. Miiller^^, 

R. Muresan28'38, B. Muryn^^, B. Muster^s, M. Musy^^, J. Mylroie-Smith^^, P. Naik^^^ 

T. Nakada^^, R. Nandakumar^^ , L Nasteva\ M. Nedos^, M. Needham^'^, N. Neufeld^^, 

C. Nguyen-Mau^^'", M. NicoF, V. Niess^, N. Nikitin3\ A. Nomerotski^^ A. Novoselov^^, 

A. Oblakowska-Mucha^^ , V. Obraztsov^'' , S. Oggero^^, S. Ogilvy'^'', O. Okhrimenko^^ , 

R. Oldeman^^''', M. Orlandea^^ J.M. Otalora Goicochea^, P. Owen^*^, K. Pal^^ J. Palacios^^ 

A. Palano^^'^, M. Palutan^^, J. Panman^'', A. Papanestis^^ , M. Pappagallo^^ , C. Parkes^^'^^, 

C.J. Parkinson^^ G. Passaleva^^ G.D. Patel^^ M. Patel^^ S.K. Paterson4^ G.N. Patrick^^ 

C. Patrignani^^'*, C. Pavel-Nicorescu^^ , A. Pazos Alvarez^^, A. Pellegrino^^ , G. Penso^^'', 
M. Pepe Altarelli^'^, S. Perazzini^^''^, D.L. Perego^Oj, E. Perez Trigo^*^, 

A. Perez-Calero Yzquierdo^^, P. Perret^, M. Perrin-Terrin^, G. Pessina^'^, A. Petrella^^'^^, 

A. Petrolini^^'S A. Phan^^^ ^ Picatoste 011oqui3^ B. Pie Valls2^ B. Pietrzyk^, T. Pilaf^^, 

D. Pinci22, R. Plackett^^ S. Playfer^'S, M. Plo Casasus^^, G. Polok^s, A. Poluektov^^'^s, 

E. Polycarpo^, D. Popov^^, B. Popovici^^, C. Potterat^^, A. Powell^\ T. du Pree^^, 
J. Prisciandaro^*^ , V. Pugatch^^, A. Puig Navarro'^^, W. Qian^^, J.H. Rademacker^^ , 

B. Rakotomiaramanana^^ , M.S. Ranged, L Raniuk^°, G. Raven^^, S. Redford^\ M.M. Reid^^, 

A. C. dos Reis\ S. Ricciardi''^ , K. Rinnert^*^, D.A. Roa Romero^, P. Robbe'^, E. Rodrigues^^'^°, 

F. Rodrigues^, P. Rodriguez Perez^^, G.J. Rogers^^, S. Roiser^'^, V. Roman ovsky^^ , 

M. Rosello^^'", J. Rouvinet3^ T. RuF, H. Ruizes, G. Sabatino^^'^ J.J. Saborido Silva^^, 
N. Sagidova^s, P. Sail'^^, B. Saitta^^''', C. Salzmann^^, M. Sannino^^-*, R. Santacesaria^^, 

C. Santamarina Rios^^, R. Santinelli^''', E. Santovetti^^''^, M. Sapunov^, A. Sarti^^'', 
C. Satriano22."^, A. Satta^^, M. Savrie^*5'% D. Savrina^o, P. Schaack'^^ M. Schiller^^, 

S. Schleich^, M. Schlupp^, M. Schmelling^°, B. Schmidt^^, O. Schneider^^, A. Schopper^^, 
M.-H. Schune^, R. Schwemmer^^, B. Sciascia^*^, A. Sciubba^^'', M. Seco^^, A. Semennikov^'^, 
K. Senderowska^^ , L Sepp^^, N. Serra^^, J. Serrano^, P. Seyfert^^, B. Shao^, M. Shapkin^^, 
L Shapoval^O'37, R Shatalov^o, Y. Shcheglov^^, T. Shears^^ L. Shekhtman^^, O. Shevchenko^° , 
V. Shevchenko^o, A. Shires''^ R. Silva Coutinho^^ T. Skwarnicki^^ A.C. Smith^'^, 
N.A. Smith'^^ E. Smith^i'^s^ k. Sobczak^, F.J.P. Soler'^^ A. Solomin^^ F. Soomrol^ 

B. Souza De Paula^, B. Spaan^, A. Sparkes^^ P. Spradlin^^ F. Stagni^^, S. Stahll^ 

O. Steinkamp^^, S. Stoica^^ S. Stone^^.s?^ ^ Storaci2^ M. Straticiuc^*^ , U. Straumann^^, 



iv 



V.K. Subbiah'^^, S. Swientek^, M. Szczekowski^'', P. Szczypka^*, T. Szumlak^^, S. T'Jampens^, 
E. Teodorescu^*^, F. Teubert^^, C. Thomas^\ E. Thomas^"^, J. van Tilburg^"*^, V. Tisserand^, 
M. Tobin^s, S. Topp-Joergensen5\ N. Torr5\ E. Tournefier'''''^, M.T. Tran^^, 
A. Tsaregorodtsev^, N. Tuning^'^, M. Ubeda Garcia^''', A. Ukleja^'', P. Urquijo^^, U. Uwer^^, 
V. Vagnoni^'*, G. Valenti^^, R. Vazquez Gomez^^, P. Vazquez Regueiro^^, S. Vecchi'^^, 
J.J. Velthuis''^, M. Veltri^^'S, B. Viaud^ I. Videau^ X. Vilasis-Cardona^s.", J. Visniakov^^, 
A. Vollhardt^^, D. Volyanskyy^^ , D. Voong^^, A. Vorobyev^^, H. Voss^'', S. Wandernoth^^, 
J. Wang52^ Ward^^ N.K. Watson5^ A.D. Webber^o, D. Websdale''^ M. Whitehead^^ 
D. Wiednerl^ L. Wiggers^^, G. Wilkinson^i, M.P. Williams^'^'^^ M. Williams''^ F.F. Wilson^^ 
J. Wishahi^ M. Witek^^, W. Witzeling^^, S.A. WotW^ K. Wyllie^^ Y. Xie^^ F. Xmg5\ 
Z. Xing^^ Z. Yang3, R. Young^^, O. Yushchenko^^, M. Zavertyaev^°''' , F. Zhang^, L. Zhang^^ 
W.C. Zhang-*^^, Y. Zhang^, A. Zhelezov^^, L. Zhong'^, E. Zverev^^, A. Zvyagin^'^. 

^ Centra Brasileiro de Pesquisas Fisicas ( CBPF), Rio de Janeiro, Brazil 
Universidade Federal do Rio de Janeiro (UFRJ), Rio de Janeiro, Brazil 
^ Center for High Energy Physics, Tsinghua University, Beijing, China 
^LAPP, Universite de Savoie, CNRS/IN2P3, Annecy-Le-Vieux, France 

^Clermont Universite, Universite Blaise Pascal, CNRS/IN2P3, LPC, Clermont-Ferrand, France 
^CPPM, Aix-Marseille Universite, CNRS/IN2P3, Marseille, France 
"^LAL, Universite Paris-Sud, CNRS/IN2P3, Orsay, France 

^ LPNHE, Universite Pierre et Marie Curie, Universite Paris Diderot, CNRS/IN2P3, Paris, France 
^ Fakultdt Physik, Technische Universitdt Dortmund, Dortmund, Germany 

Max- Planck- Institut fiir Kernphysik (MPIK), Heidelberg, Cermany 

Physikalisches Institut, Ruprecht-Karls- Universitdt Heidelberg, Heidelberg, Germany 

School of Physics, University College Dublin, Dublin, Ireland 
^^Sezione INFN di Bari, Bari, Italy 
^'^Sezione INFN di Bologna, Bologna, Italy 

Sezione INFN di Cagliari, Cagliari, Italy 

Sezione INFN di Ferrara, Ferrara, Italy 

Sezione INFN di Firenze, Firenze, Italy 

Laboratori Nazionali dell'INFN di Frascati, Frascati, Italy 

Sezione INFN di Cenova, Genova, Italy 
'^^Sezione INFN di Milano Bicocca, Milano, Italy 
^^Sezione INFN di Roma Tor Vergata, Roma, Italy 
'^'^ Sezione INFN di Roma La Sapienza, Roma, Italy 

'^^Nikhef National Institute for Subatomic Physics, Amsterdam, The Netherlands 

'^'^Nikhef National Institute for Subatomic Physics and Vrije Universitdt, Amsterdam, The Netherlands 
^^Henryk Niewodniczanski Institute of Nuclear Physics Polish Academy of Sciences, Kracow, Poland 
^^AGH University of Science and Technology, Kracow, Poland 

Soltan Institute for Nuclear Studies, Warsaw, Poland 
^^Horia Hulubei National Institute of Physics and Nuclear Engineering, Bucharest-Magurele, Romania 

Petersburg Nuclear Physics Institute (PNPI), Catchina, Russia 

Institute of Theoretical and Experimental Physics (ITEP), Moscow, Russia 

Institute of Nuclear Physics, Moscow State University (SINP MSU), Moscow, Russia 
'^'^ Institute for Nuclear Research of the Russian Academy of Sciences (INR RAN), Moscow, Russia 
^^Budker Institute of Nuclear Physics (SB RAS) and Novosibirsk State University, Novosibirsk, Russia 
'^'^ Institute for High Energy Physics (IHEP), Protvino, Russia 

Universitat de Barcelona, Barcelona, Spain 

Universidad de Santiago de Compostela, Santiago de Compostela, Spain 
'^'^ European Organization for Nuclear Research ( CERN), Geneva, Switzerland 
'^^Ecole Polytechnique Federale de Lausanne (EPFL), Lausanne, Switzerland 



V 



Physik-Institut, Universitdt Zurich, Zurich, Switzerland 
'^^ NSC Kharkiv Institute of Physics and Technology (NSC KIPT), Kharkiv, Ukraine 

Institute for Nuclear Research of the National Academy of Sciences (KINR), Kyiv, Ukraine 
'^^H.H. Wills Physics Laboratory, University of Bristol, Bristol, United Kingdom 

Cavendish Laboratory, University of Cambridge, Cambridge, United Kingdom 
^'^Department of Physics, University of Warwick, Coventry, United Kingdom 

STFC Rutherford Appleton Laboratory, Didcot, United Kingdom 

School of Physics and Astronomy, University of Edinburgh, Edinburgh, United Kingdom 
^"^ School of Physics and Astronomy, University of Clasgow, Glasgow, United Kingdom 

Oliver Lodge Laboratory, University of Liverpool, Liverpool, United Kingdom 
'^^ Imperial College London, London, United Kingdom 

^'^ School of Physics and Astronomy, University of Manchester, Manchester, United Kingdom 
Department of Physics, University of Oxford, Oxford, United Kingdom 
Syracuse University, Syracuse, NY, United States 

CC-IN2P3, CNRS/IN2P3, Lyon- Villeurbanne, France, associated member 
^'^Pontificia Universidade Catolica do Rio de Janeiro (PUC-Rio), Rio de Janeiro, Brazil, associated to ^ 
University of Birmingham, Birmingham, United Kingdom 

"^P.N. Lebedev Physical Institute, Russian Academy of Science (LPI RAS), Moscow, Russia 



^ Universitd di 
Universitd di 
Universitd di 
Universitd di 



* Universitd di 
^ Universitd di 
Universitd di 



Bari, Bari, Italy 
Bologna, Bologna, Italy 
Cagliari, Cagliari, Italy 
Ferrara, Ferrara, Italy 



^ Universitd di Firenze, Firenze, Italy 

3 Universitd di Urbino, Urbino, Italy 

^ Universitd di Modena e Reggio Emilia, Modena, Italy 



Cenova, Genova, Italy 
Milano Bicocca, Milano, Italy 
Roma Tor Vergata, Roma, Italy 



' Universitd di Roma La Sapienza, Roma, Italy 
™ Universitd della Basilicata, Potenza, Italy 

^LIFAELS, La Salle, Universitat Ramon Llull, Barcelona, Spain 
"Hanoi University of Science, Hanoi, Viet Nam 



vi 



1 Introduction 



The charged particle multiphcity is a basic observable that characterizes the hadronic final 
state. The multiplicity distribution is sensitive to the underlying QCD dynamics of the 
proton-proton collision. ALICE [ij, ATLAS jij and CMS jsj have measured the charged 
multiplicity distributions mainly covering the central region, while LHCb's geometrical 
acceptance allows the dynamics of the collision to be probed in the forward region. The 
forward region is in particular sensitive to low Bjorken-x QCD dynamics and multi-parton 
interactions (MPI) i]. 

In this analysis, the charged particles are reconstructed in the vertex detector (VELO) 
surrounding the interaction region. The VELO was designed to provide a uniform accep- 
tance in the forward region with additional coverage of the backward region. In the 
absence of almost any magnetic field in the VELO region, the particle trajectories are 
straight lines and therefore no acceptance corrections as a function of momentum are 
needed. Since the VELO is close to the interaction region, the amount of material before 
the particle detection is small, minimising the corrections for particle interactions with 
detector material. 

This paper is organized as follows. Section [2] gives a brief description of the LHCb 
detector and the configuration used to record data in Spring 2010. The Monte Carlo 
simulation and data selection are outlined in Sections [3] and H] respectively, with Sec- 
tion [S] giving an overview of the analysis. The systematic uncertainties are outlined in 
Section [61 The final results are discussed in Section [7] and compared with different model 
expectations, before concluding in Section [HI 



2 LHCb detector 

The LHCb detector is a single-arm magnetic dipole spectrometer with a polar angular 
coverage with respect to the beam line of approximately 15 to 300 mrad in the horizontal 
bending plane, and 15 to 250 mrad in the vertical non-bending plane. The detector is 
described in detail elsewhere [l]. A right-handed coordinate system is defined with its 
origin at the nominal proton-proton interaction point, the z axis along the beam line and 
pointing towards the magnet, and the y axis pointing upwards. 

For the low luminosity running period of the LHC relevant for this analysis, the prob- 
ability of observing more than one collision in a proton-proton bunch crossing (pile-up) is 
measured to be (3.7 ± 0.4)%, dominated by a double interaction. For the measurements 
presented in this paper the tracking detectors are of particular importance. The LHCb 
tracking system consists of the VELO surrounding the proton-proton interaction region, a 
tracking station (TT) before the dipole magnet, and three tracking stations (T1-T3) after 
the magnet. Particles traversing from the interaction region to the downstream tracking 
stations experience an integrated bending-field of approximately 4 Tm. 

The VELO consists of silicon microstrip modules, providing a measure of the radial 
and azimuthal coordinates, r and 0, distributed in 23 stations arranged along the beam 
direction. The first two stations at the most upstream z positions are instrumented to 
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provide information on the number of visible interactions in the detector at the first level 
of the trigger. The VELO is constructed in two halves, movable in the x and y directions 
so that it can be centered on the beam. During stable beam conditions the two halves 
are located at their nominal closed position, with active silicon at only 8 mm from the 
beams, providing full azimuthal coverage. 

The TT station also uses silicon microstrip technology. The T1-T3 tracking stations 
have silicon microstrips in the region close to the beam pipe, whereas straw tubes are 
employed in the outer region. 

Though the particle multiplicity is measured using only tracks reconstructed with the 
VELO, momentum information is only available for "long" tracks. Long tracks are formed 
from hits in the VELO (before the magnet) and in the T1-T3 stations (after the magnet). 
If available, measurements in the TT station are added to the long track. 

The LHCb trigger system consists of two levels. The first level is implemented in 
hardware and is designed to reduce the event rate to a maximum of 1 MHz. The complete 
detector is then read out and the data is sent to the second level, a software trigger. 
For the early data taking period with low luminosity used in this analysis a simplified 
trigger was used. The first level trigger was operated in pass-through mode. A fast track 
reconstruction was performed in the software trigger and events with at least one track 
observed in the VELO were accepted. 



3 Monte Carlo simulation 

Monte Carlo event simulation is used to correct for acceptance, resolution effects and 
for background characterisation. The detector simulation is based on the Geant4 
package. Details of the detector simulation are given in Ref. [5|. The simulated material 
of the components of the VELO was compared with the masses measured at the time of 
production and agreement was found to be within 15%. The Monte Carlo event samples 
are passed through reconstruction and selection procedures identical to those for the data. 

Elastic and inelastic proton-proton collisions are generated using the Pythia 6.4 event 
generator [3], with CTEQ6L parton density functions j^, which is tuned to lower energy 
hadron collider data j^. The inelastic processes include both single and double diffractive 



components. The decay of the generated particles is carried out by EvtGen [10|, with 



final state radiation handled by Photos Secondary particles produced in material 



interactions are decayed through the Geant4 program. 



4 Data selection 

A sample of 3 x 10^ events, collected during May 2010, was used in this analysis. In order 
to minimize the contribution of secondary particles and misreconstructed (fake) tracks, 
only the tracks satisfying a set of minimal quality criteria are accepted. To minimise fake 
tracks a cut on the P^^' degree of freedom of the reconstructed track, x^/ndf < 5, 
is applied. To further reduce fake tracks, and reduce duplicate tracks due to a split of 
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the reconstructed trajectory, a cut of less than four missing VELO hits compared to the 
expectation is apphed. To ensure that tracks originate from the primary interaction, the 
requirements do < 2 mm and zq < Sai are apphed, where do is the track's closest distance 
to the beam line, zq is the distance along the z direction from the centre of the luminous 
region and ai is the width of the luminous region extracted from a Gaussian fit. 

Tracks are considered for this analysis only if their pseudorapidity is in either of the 
ranges —2.5 < r] < —2.0 and 2.0 < r] < 4.5. Pseudorapidity is defined as — ln[tan(6'/2)] 
where 6 is the polar angle of the particle with respect to the z direction. The forward 
range is divided in five equal sub-intervals with Arj = 0.5. 

5 Analysis strategy 

The reconstructed multiplicity distributions are corrected on an event by event basis to 
account for the tracking and selection efficiencies and for the background contributions. 
These corrected distributions are then used to measure the charged particle multiplici- 
ties in each of the t] intervals (bins) through an unfolding procedure. Only events with 
tracks in the rj bins are included in the distributions and subsequent normalisation. The 
distributions are corrected for pile-up effects so they represent the charged particle multi- 
plicities, Uch, for single proton-proton interactions. No unfolding procedure is required for 
the charged particles pseudorapidity density distribution i.e. the mean number of charged 
particles per single pp-collision and unit of pseudorapidity. Only the per track corrections 
for background and tracking efficiency are needed. For this distribution, at least one 
VELO track is required in the full forward rj range. Each of these elements of the analysis 
procedure are discussed in subsequent subsections. 

Hard interaction events are defined by requiring at least one long track with px > 
1 GeV/c in the range 2.5 < rj < 4.5 where the detector has high efficiency. The geometric 
acceptance is no longer independent of momentum and therefore the distributions require 
an additional correction. 

In this analysis primary charged particles are defined as all particles for which the 
sum of the ancestors' mean lifetimes is smaller than 10 ps; according to this definition the 
decay products of beauty and charm are primary particles. 

5.1 Efficiency correction 

The LHCb simulation is used to estimate the overall tracking and selection efficiency as a 
function of pseudorapidity and azimuthal angle (p. As the VELO is outside the magnetic 
field region tracks are straight lines and no study of acceptance as a function of momentum 
is necessary. It is found that the efficiency (including acceptance) in the forward region 
is typically greater than 90% while it is at least 85% in the backward region. Tracking 
efficiency depends weakly on the event track multiplicity; this is taken into account in the 
evaluation of the systematic error. 
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5.2 Background contributions 



There are two main sources of background that can affect the measurement of the muhi- 
phcity of charged particles: secondary particles misidentified as primary and fake tracks. 
Other sources of background, such as beam-gas interactions, are estimated to be negligi- 
ble. 

The correlation between the number of VELO hit clusters in an event and its track 
multiplicity is in good agreement between the data and simulation, indicating that the 
fraction of fake tracks is well understood. It is also found that for each rj bin the multi- 
plicity of fake tracks is linearly dependent on the number of VELO clusters in the event. 
Therefore it is possible to parameterise the fake contribution as a function of VELO 
clusters using the Monte Carlo simulation. 

The majority of secondary particles are produced in photon conversions in the VELO 
material, and in the decay of long-lived strange particles such as Kg and hyperons. While 
earlier LHCb measurements show that the production of is reasonably described by 



the Monte Carlo generator [12|, there are indications that the production of A particles 
is underestimated 13|] . This difference is accounted for in the systematic error associated 
with the definition of primary particles. 

The fraction of secondary particles is estimated as a function of both rj and (p. In 
general, depending on the t] bin, the correction for non-primary particles (from conversion 
and secondaries) changes the mean values of the particle multiplicity distributions by 
5 - 10%. 



5.3 Correction and unfolding procedure 

The procedure consists of three steps; a background subtraction is made, followed by an 
efficiency correction and finally a correction for pile-up. The procedure is applied to all 
measured track multiplicity distributions in each of the different rj intervals. 

In the first step, the distribution is corrected for fake tracks and non-primary particles. 
A mean number of background tracks is estimated for each event based on the parame- 
terizations described in Section 15.21 A PDF (probability density function) is built with 
this mean value assuming a Poisson distribution for the number of background tracks. 
Hence, a PDF for the number of prompt charged particles in a given event is then ob- 
tained. These per event PDFs are summed up and normalized to obtain the reconstructed 
prompt charged track multiplicity distribution i. e. the fraction of events with utr tracks, 
Prob(nti.). 

In the second step, the correction for the tracking efficiency is applied. For each rj bin 
a mean efficiency, e, is calculated based on the per track efficiency as function of (^7,0). 
As explained below, this is used to unfold the background-subtracted track multiplicity 
distribution, Prob(nti.), to obtain the underlying charged particle multiplicity distribution, 
Prob(n(;h), where rich is the number of primary produced particles of all proton-proton 
collisions in an event. 

For a given value of fich, the probability to observe Utr reconstructed tracks given a 
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reconstruction efficiency e is described by the binomial distribution 

P(ntn iich, e) = ( ) (1 - e)^--°-6°-. (1) 
Hence, the observed track multiphcity distribution is given by 

oo 

Prob(ntr) = Prob(nch) x p(ntr, n<.h, e). (2) 

nch=0 

The values for Prob(nch) are obtained by performing a fit to Prob(ntr). The procedure 
has been verified using simulated data. 

In the last step, the distributions are corrected for pile-up to obtain the charged par- 
ticle multiplicity distributions of single interaction events, Prob(nch). This is done using 
an iterative procedure. For low luminosity, Prob(iich) has mainly two contributions: sin- 
gle proton-proton interactions and a convolution of two single proton-proton interactions. 
The starting assumption is that the observed distribution is the single proton-proton in- 
teraction. From this, the convolution term is calculated, and by subtracting it from the 
observed distribution, a ffist order estimate for the single proton-proton distribution is ob- 
tained. This can then be used to calculate again the convolution term and obtain a second 
order estimate for the single proton-proton distribution. The procedure usually converges 
after the second iteration. The pile-up correction typically changes the mean value of the 
particle multiplicity distributions by 3 — 4%. It was checked that the contribution from 
pile-up events with more than two proton-proton collisions is negligible. 

As mentioned before, no unfolding procedure is required for the charged particles 
pseudorapidity density, only the per track corrections for background tracks and track- 
ing efficiency are applied. The distribution is then normalized to the total number of 
proton-proton collisions including pile-up collisions. In the case of hard interactions, the 
pseudorapidity density distribution of the pile-up collisions without the cut is ffist 
subtracted. Finally, the distribution is normalized to the total number of hard collisions. 



6 Systematic uncertainty 
6.1 Efficiency 

Studies based on data and simulation show that the error on the tracking efficiency for 
particles reaching the tracking stations T1-T3 is < 3% The tracking efficiency re- 
duces for low- momentum (px < 50MeV/c) particles due to interactions with the detector 
material and the residual magnetic field in the VELO region. Since no momentum mea- 
surement exists for the reconstructed VELO tracks, the estimate of a mean efficiency relies 
on the prediction of the LHCb Monte Carlo model for the contribution of low-momentum 
particles to the total number of particles. The simulation predicts that in the forward 
region the fraction of particles below a transverse momentum of 50MeV/c is 2.4%. The 
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corresponding average single track efficiency in this rj range is measured to be 94%. In 
the two extreme cases in which no particles with below 50MeV/c were reconstructed 
or no such particles were produced the average track efficiency would be reduced by 1.2% 
or increased by 1.1% respectively. Assuming a 25% uncertainty on the number of low 
momentum particles, as suggested by the comparison between the measured particle mul- 
tiplicity and Monte Carlo prediction, the additional contribution to the track efficiency 
uncertainty is < 1%. Adding this to the 3% track reconstruction uncertainty, gives an 
overall 4% error on the track efficiency used in the unfolding procedure. The systematic 
error contribution is then estimated by unfolding the multiplicity distributions varying 
the tracking efficiency by ±4%. 

6.2 Non-primary particles 

The main systematic uncertainty on the contribution of non-primary particles arises from 
the knowledge of the detector material (15%). Two thirds of non-primary particles are 
due to conversions of photons from vr'' decays, resulting in an 10% uncertainty. The 
multiplicity of vr*^ scales with the charged multiplicity, therefore no additional error is 
applied. Varying by ±40% the production of A results in an uncertainty of about 5% on 
the non-primary contribution. A pessimistic assumption of a 25% underestimation of the 
non-prompt contribution would change the mean and RMS values of the particle multi- 
plicity distributions by —2%, which can be neglected compared to the tracking efficiency 
uncertainty of 4%. 

6.3 Pile-up 

The pile-up corrections inherit a systematic uncertainty from the determination of the 
mean number of visible interactions of 10%. This correction to the pile-up fraction is 
small and is negligible compared to the systematic uncertainty due to the track efficiency 
correction. 



7 Results 

Figure [T] shows the unfolded charged particle multiplicity distribution for different bins 
in pseudorapidity, t]. Figure [2] shows the multiplicity distributions for the full forward 
range, 2.0 < rj < 4.5. There is a requirement of at least one track in the relevant t] range. 
The distributions are compared to several Monte Carlo event generators. Pythia 6.424 
is compared with the data for a number of tunes including the LHCb tuned settings joj. 



In particular the PerugiaO and PerugiaNOCR tunings [15[ are shown. In addition, the 



Pythia 8.145 generator [16| was compared to the data as well as Phojet 1-12.35 [17[. In 
general all generators underestimate the multiplicity distributions, with the LHCb tune 
giving the best description of the data; this tune does not use data from the LHC. The 
exclusion of the Pythia diffractive processes in the Perugia tunes. Figs. [T|d and [2]3, also 
improves the description of the data, particularly in the full forward region. 
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"ch "ch 

Figure 1: The multiplicity distribution in r] bins (shown as points with statistical error bars) 
with predictions of different event generators. The inner error bar represents the statistical 
uncertainty and the outer error bar represents the systematic and statistical uncertainty on the 
measurements. The data in both figures are identical with predictions from Pythia 6, Phojet 
and Pythia 8 in (a) and predictions of the Pythia 6 Perugia tunes with and without diffraction 
in (b). 



The Koba-Nielsen-Olesen (KNO) scaling variable [18| has been used to compare the 
data in the different t] bins. Figure E] shows the KNO scaled multiplicity distributions, 
"^{u) = (uch) X Prob(nch) as a function of m = As the multiplicity distributions 

measured are truncated the mean used was extracted by fitting a negative binomial dis- 
tribution. It clearly shows that the distributions in the different rj bins are equivalent. In 
particular this illustrates that when there is a requirement of at least one track in the t] 
bin the forward and backward regions (2.0 < |?7| < 2.5) are identical. 

The charged particle pseudorapidity density, p, is shown as a function of pseudorapidity 
in Fig. m The data have a marked asymmetry between the forward and backward region; 
this is a consequence of the requirement of at least one track in the full forward rj range. All 
models fail to describe the mean charged particle multiplicity per unit of pseudorapidity. 
The models, to varying degrees, also display the asymmetry but in none of the models 
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Figure 2: The multiplicity distribution in the forward i] range (shown as points with error bars) 
with predictions of different event generators. The shaded bands represent the total uncertainty 
on the measurements. The data in both figures are identical with predictions from Pythia 6, 
Phojet and Pythia 8 in (a) and predictions of the Pythia 6 Perugia tunes with and without 
diffraction in (b). 

is this as large as in the data. The effect on the predictions of excluding diffractive 
processes is shown in Fig. |4)3 using the Perugia tunes. There is a better description of the 
rj distribution in the backward directions but it still fails to describe the forward-backward 
asymmetry. 

A sample of hard QCD events were studied by ensuring at least one track in the 
pseudorapidity range 2.5 < t] < 4.5 has a transverse momentum > IGeV/c. In 
comparison to the data without this pt requirement, the multiplicity distributions have 
larger high multiplicity tails, see Figs. |5]and|6l The data are again compared to predictions 
of several event generators. In general the predictions are in better agreement than for the 
minimum bias data but the pseudorapidity range 4.0 < t] < 4.5 remains poorly described. 
As the Pt cut removes the majority of diffractive events from Pythia 6 the comparisons 
with and without diffraction are not shown. 
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Figure 3: The KNO distributions in different bins of rj. Only the the statistical uncertainties 
are shown. 

The charged particle density as a function of pseudorapidity for the hard QCD sample 
is shown in Fig. [71 The discontinuity observed in the data at 77 = 2.5 is an artefact 
of the event selection for the hard events. The asymmetry between the forward and 
backward region is further amplified in this sample. All models fail to describe the mean 
charged particle multiplicity per unit of pseudorapidity. The models, to varying degrees, 
also display the asymmetry but never give an effect as large as the data. The Perugia 
(NOCR) tune gives the best description of the data in the backward direction but fails 
to reproduce the size of the asymmetry. 

8 Summary 

The LHCb spectrometer acceptance, 2.0 < rj < 4.5, allows the forward region to be 
probed at the LHC. The charged multiplicity distributions at ^/s = 7 TeV are measured 
with and without a pt event selection, making use of the high efficiency of the LHCb 
VELO. Several event generators are compared to the data; none are fully able to describe 
the multiplicity distributions or the charged density distribution as a function of 77 in the 
LHCb acceptance. In general, the models underestimate the charged particle production, 
in agreement with the measurements in the central region at the LHC. 
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Figure 4: The charged particle densities as a function of r/ (shown as points with statistical error 
bars) and comparisons with predictions of event generators, as indicated in the key. The shaded 
bands represent the total uncertainty. The events are selected by requiring at least one charged 
particle in the range 2.0 < r] < 4.5. The data in both figures are identical with predictions from 
Pythia 6, Phojet and Pythia 8 in (a) and predictions of the Pythia 6 Perugia tunes with 
and without diffraction in (b). 
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Figure 5: The multiplicity distribution in rj bins (shown as points with error bars) with pre- 
dictions of different event generators. The inner error bar represents the statistical uncertainty 
and the outer error bar represents the systematic and statistical uncertainty on the measure- 
ments. The events have at least one track with a pT > l.OGeV/c in the pseudorapidity range 
2.5 < rj < 4.5. The data in both figures are identical with predictions from Pythia 6, Phojet 
and Pythia 8 in (a) and predictions of the Pythia 6 Perugia tunes in (b). 
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Figure 6: The multiplicity distribution in the forward rj range (shown as points with statistical 
error bars) with predictions of different event generators. The shaded bands represent the total 
uncertainty. The events have at least one track with a pi > l.OGeV/c in the pseudorapidity 
range 2.5 < i] < 4.5. The data in both figures are identical with predictions from Pythia 6, 
Phojet and Pythia 8 in (a) and predictions of the Pythia 6 Perugia tunes in (b). 
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Figure 7: The data charged particle densities as a function of rj (shown as points with statistical 
error bars) and comparisons with predictions of event generators, as indicated in the key. The 
events have at least one track with a pT > l.OGeV/c in the pseudorapidity range 2.5 < r] < 4.5. 
The shaded bands represent the total uncertainty. 
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A Tables of charged particle multiplicities 



Table 1: Charged particle multiplicity distribution in the pseudorapidity range —2.5 < rj < —2.0 
for minimum bias events and for hard QCD events (see text). The first quoted uncertainty is 
statistical and the second is systematic. 
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60 83 ± 45 ± 1 13 


79 84 ± 79 ± 48 


7 


46.08 ± 0.48 ± 1.33 


63.40 ± 0.83 ± 1.33 


8 


35.01 ± 0.50 ± 1.35 


51.30 ± 0.90 ± 1.63 


9 


26.43 ± 0.52 ± 1.40 


40.66 ± 0.97 ± 1.81 


10 


19.75 ±0.55 ±1.36 


31.50 ± 1.02 ± 1.86 


11 


14.60 ± 0.57 ± 1.19 


24.16 ± 1.08 ± 1.83 


12 


10.82 ±0.59 ±1.00 


18.03 ± 1.12 ± 1.64 


13 


7.86 ±0.61 ±0.90 


13.96 ± 1.21 ± 1.61 


14 


5.57 ±0.63 ±0.86 


9.56 ± 1.19 ± 1.28 


15 


3.94 ±0.65 ±0.73 


7.14 ± 1.30 ± 1.09 


16 


2.90 ±0.67 ±0.37 


5.10 ± 1.29 ± 1.11 


17 


2.44 ± 0.68 ± 0.96 


4.48 ± 1.34 ± 1.28 


18 


1.14 ±0.70 ±0.61 


2.13 ± 1.43 ±2.03 


19 


0.96 ±0.71 ±0.66 


1.78 ±1.41 ±0.19 


20 


0.75 ±0.72 ±0.27 


1.46 ± 1.44 ±0.60 
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Table 2: Charged particle multiplicity distribution in the pseudorapidity range 2.0 < r] < 2.5 
for minimum bias events and for hard QCD events (see text). The first quoted uncertainty is 
statistical and the second is systematic. 





Prob. in min. bias 


Prob. in hard QCD 




events xlO^ 


events x 10'^ 


1 


244.35 ± 0.36 ± 7.66 


126.88 ±0.38 ±6.57 


2 


191.00 ±0.33 ±4.02 


140.50 ±0.43 ±5.81 


3 


142.72 ±0.31 ±1.44 


133.83 ±0.44 ±3.91 


4 


106.75 ±0.28 ±0.10 


121.45 ±0.44 ±1.95 


5 


80.27 ±0.26 ±0.73 


103.10 ±0.43 ±0.75 


6 


61.09 ± 0.25 ± 1.22 


86.87 ±0.42 ±0.98 


7 


46.22 ±0.23 ±1.42 


70.01 ± 0.41 ± 1.59 


8 


34.57 ±0.21 ± 1.45 


55.15 ± 0.39 ± 1.83 


9 


26.09 ±0.20 ±1.38 


43.12 ±0.36 ±2.13 


10 


19.30 ±0.18 ±1.34 


32.71 ±0.34 ±2.20 


11 


14.08 ± 0.17 ± 1.17 


24.64 ±0.32 ±2.00 


12 


10.17± 0.16 ±1.07 


18.25 ± 0.29 ± 1.80 


13 


7.23 ±0.14 ±0.98 


13.66 ± 0.28 ± 1.84 


14 


5.43 ±0.13 ±0.82 


9.97 ± 0.25 ± 1.52 


15 


3.55 ±0.12 ±0.60 


6.64 ± 0.22 ± 1.12 


16 


2.60 ±0.11 ±0.40 


4.91 ±0.21 ±0.78 


17 


1.78 ±0.10 ±0.65 


3.14 ± 0.18 ± 1.23 


18 


1.35 ±0.09 ±0.28 


2.45 ±0.17 ±0.47 


19 


0.82 ±0.08 ±0.22 


1.56 ±0.15 ±0.42 


20 


0.62 ±0.07 ±0.19 


1.15 ±0.13 ±0.34 
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Table 3: Charged particle multiplicity distribution in the pseudorapidity range 2.5 < r] < 3.0 
for minimum bias events and for hard QCD events (see text). The first quoted uncertainty is 
statistical and the second is systematic. 



Thrh 


Prob. in min. bias 


Prob. in hard QCD 




events xlO^ 


events x 10'^ 


1 


249.37 ±0.35 ±7.88 


121.02 ±0.36 ±6.72 


2 


194.45 ±0.33 ±4.11 


140.71 ±0.41 ±6.20 


3 


144.53 ± 0.29 ± 1.39 


138.90 ±0.42 ±4.26 


4 


107.18 ±0.27 ±0.10 


125.71 ±0.41 ±2.10 


5 


80.42 ± 0.24 ± 0.89 


108.13 ±0.40 ±0.34 


6 


60.29 ± 0.22 ± 1.34 


87.75 ± 0.37 ± 1.24 


7 


45.03 ± 0.20 ± 1.53 


70.69 ± 0.35 ± 1.85 


8 


33.53 ± 0.18 ± 1.55 


55.79 ±0.33 ±2.31 


9 


24.75 ±0.16 ±1.46 


42.12 ±0.30 ±2.40 


10 


17.98 ±0.15 ±1.30 


31.82 ±0.27 ±2.23 


11 


12.98 ± 0.13 ± 1.23 


23.37 ±0.25 ±2.10 


12 


9.16 ± 0.12 ± 1.12 


16.64 ± 0.22 ± 1.95 


13 


6.74 ±0.11 ±0.87 


12.07 ± 0.19 ± 1.52 


14 


4.46 ±0.09 ±0.71 


8.43 ±0.17 ±1.27 


15 


3.23 ±0.08 ±0.47 


5.97 ±0.15 ±0.88 


16 


2.20 ±0.07 ±0.71 


4.07 ± 0.13 ± 1.31 


17 


1.57 ±0.06 ±0.32 


2.78 ±0.11 ±0.52 


18 


0.94 ±0.05 ±0.32 


1.86 ±0.10 ±0.51 


19 


0.69 ±0.05 ±0.33 


1.26 ±0.09 ±0.56 


20 


0.50 ±0.04 ±0.13 


0.92 ±0.08 ±0.20 
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Table 4: Charged particle multiplicity distribution in the pseudorapidity range 3.0 < r] < 3.5 
for minimum bias events and for hard QCD events (see text). The first quoted uncertainty is 
statistical and the second is systematic. 



Thrh 


Prob. in min. bias 


Prob. in hard QCD 




events xlO^ 


events x 10'^ 


1 


257.54 ±0.36 ±8.38 


128.89 ±0.38 ±7.33 


2 


199.12 ±0.33 ±4.08 


145.79 ±0.41 ±6.39 


3 


147.50 ± 0.30 ± 1.23 


145.41 ±0.43 ±4.13 


4 


108.21 ±0.27 ±0.31 


130.01 ±0.42 ±2.16 


5 


79.83 ± 0.24 ± 1.10 


109.73 ±0.41 ±0.44 


6 


58.83 ± 0.22 ± 1.50 


87.48 ± 0.38 ± 1.58 


7 


43.25 ± 0.20 ± 1.67 


67.91 ±0.35 ±2.16 


8 


31.48 ± 0.18 ± 1.64 


52.94 ±0.32 ±2.50 


9 


22.72 ±0.16 ±1.48 


38.50 ±0.29 ±2.43 


10 


16.12 ± 0.14 ± 1.28 


28.21 ±0.26 ±2.21 


11 


11.37±0.13±1.19 


20.63 ±0.24 ±2.17 


12 


7.89 ±0.11 ±1.07 


14.74 ±0.21 ± 1.83 


13 


5.63 ±0.10 ±0.81 


10.02 ± 0.18 ± 1.45 


14 


3.54 ±0.08 ±0.67 


7.00 ±0.16 ±1.02 


15 


2.53 ±0.07 ±0.71 


4.49 ± 0.13 ± 1.37 


16 


1.79 ±0.06 ±0.38 


3.33 ±0.12 ±0.64 


17 


1.07 ±0.06 ±0.29 


1.96 ±0.10 ±0.53 


18 


0.75 ±0.05 ±0.17 


1.38 ±0.09 ±0.32 


19 


0.49 ±0.04 ±0.22 


0.94 ± 0.08 ± 0.43 


20 


0.35 ±0.04 ±0.10 


0.65 ±0.07 ±0.17 
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Table 5: Charged particle multiplicity distribution in the pseudorapidity range 3.5 < r] < 4.0 
for minimum bias events and for hard QCD events (see text). The first quoted uncertainty is 
statistical and the second is systematic. 





Prob. in min. bias 


Prob. in hard QCD 




events xlO^ 


events x 10'^ 


1 


268.35 ±0.37 ±8.77 


139.99 ±0.39 ±7.61 


2 


206.16 ±0.34 ±4.00 


158.42 ±0.44 ±6.72 


3 


150.62 ±0.31 ±0.98 


151.42 ±0.45 ±4.01 


4 


108.81 ±0.28 ±0.56 


133.07 ±0.44 ±1.67 


5 


78.99 ± 0.25 ± 1.35 


110.17 ±0.42 ±0.92 


6 


56.92 ± 0.22 ± 1.77 


84.74 ± 0.38 ± 1.91 


7 


40.49 ± 0.20 ± 1.81 


65.65 ±0.36 ±2.61 


8 


28.60 ± 0.18 ± 1.68 


48.06 ±0.32 ±2.71 


9 


19.98 ±0.16 ±1.46 


34.60 ± 0.29 ± 2.49 


10 


13.79 ± 0.14 ± 1.30 


24.49 ± 0.26 ± 2.26 


11 


9.31 ±0.12 ± 1.18 


16.62 ±0.22 ±2.05 


12 


6.48 ±0.11 ±0.94 


11.50 ± 0.19 ± 1.51 


13 


4.02 ±0.09 ±0.68 


7.40 ± 0.17 ± 1.18 


14 


2.80 ±0.08 ±0.41 


5.09 ±0.15 ±0.75 


15 


1.82 ±0.07 ±0.64 


3.48 ± 0.13 ± 1.27 


16 


1.24 ±0.06 ±0.28 


2.23 ±0.11 ±0.45 


17 


0.68 ±0.05 ±0.25 


1.35 ±0.09 ±0.43 


18 


0.50 ±0.04 ±0.21 


0.85 ±0.08 ±0.47 


19 


0.27 ±0.04 ±0.05 


0.55 ±0.06 ±0.14 


20 


0.18 ±0.03 ±0.08 


0.31 ±0.05 ±0.18 
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Table 6: Charged particle multiplicity distribution in the pseudorapidity range 4.0 < r] < 4.5 
for minimum bias events and for hard QCD events (see text). The first quoted uncertainty is 
statistical and the second is systematic. 





Prob. in min. bias 


Prob. in hard QCD 




events xlO^ 


events x 10'^ 


1 


284.08 ±0.40 ±9.11 


159.68 ±0.01 ±8.81 


2 


215.09 ±0.38 ±4.25 


174.85 ±0.01 ±6.65 


3 


155.18 ±0.35 ±0.72 


159.67 ±0.01 ±3.42 


4 


109.77 ± 0.32 ± 1.07 


135.15 ±0.01 ±0.61 


5 


76.74 ± 0.29 ± 1.76 


107.91 ±0.01 ± 1.45 


6 


53.34 ± 0.27 ± 1.97 


82.45 ±0.01 ±2.49 


7 


36.49 ± 0.24 ± 1.93 


58.82 ±0.01 ±2.84 


8 


24.57 ± 0.22 ± 1.75 


41.25 ±0.01 ±2.75 


9 


16.30 ±0.20 ±1.50 


28.48 ±0.01 ±2.55 


10 


10.63 ± 0.17 ± 1.25 


18.52 ±0.01 ±2.11 


11 


6.76 ± 0.15 ± 1.00 


12.41 ±0.01 ± 1.83 


12 


4.20 ±0.13 ±0.70 


7.64 ±0.01 ± 1.25 


13 


2.92 ±0.12 ±0.57 


5.63 ±0.01 ±1.12 


14 


1.48 ±0.10 ±0.86 


2.66 ±0.01 ±1.54 


15 


1.15 ±0.09 ±0.33 


2.35 ±0.01 ±0.67 


16 


0.55 ±0.07 ±0.21 


1.08 ±0.01 ±0.40 


17 


0.35 ±0.06 ±0.28 


0.71 ±0.01 ±0.54 


18 


0.24 ±0.05 ±0.12 


0.45 ±0.01 ±0.21 


19 


0.09 ±0.04 ±0.13 


0.17 ±0.01 ±0.24 


20 


0.07 ±0.04 ±0.02 


0.14 ±0.01 ±0.05 
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Table 7: Charged particle multiplicity distribution in the pseudorapidity range 2.0 < rj < 4.5 
for minimum bias events and for hard QCD events (see text). The first quoted uncertainty is 
statistical and the second is systematic. 
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8.66 ± 0.25 ± 0.62 


20.40 ± 0.63 ± 1.45 


35 


7.43 ±0.26 ±0.60 


17.59 ± 0.63 ± 1.52 


37 


6.19 ±0.26 ±0.72 


15.85 ± 0.66 ± 1.88 


39 


5.56 ±0.26 ±0.71 


13.11 ±0.64 ±1.45 


41 


4.40 ±0.25 ±0.62 


11.22 ±0.64 ±1.32 


43 


3.71 ±0.25 ±0.56 


9.55 ±0.63 ±1.24 


45 


3.14 ±0.24 ±0.44 


7.74 ± 0.59 ± 1.27 


47 


2.68 ±0.23 ±0.46 


6.21 ±0.58 ± 1.40 


49 


2.00 ±0.22 ±0.49 


5.38 ± 0.54 ± 1.09 


51 


1.70 ±0.12 ±0.32 


4.18 ± 0.30 ± 1.09 


54 


1.22 ±0.11 ±0.24 


3.04 ±0.27 ±0.69 


57 


0.88 ±0.09 ±0.20 


2.26 ±0.24 ±0.49 


60 


0.63 ±0.08 ±0.15 


1.58 ±0.21 ±0.45 
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